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From Solar and Stellar Flares to Coronal Heating: 
Theory and Observations of How Magnetic 
Reconnection Regulates Coronal Conditions 

P. A. Cassak 1 , D. J. Mullan 1 and M. A. Shay 1 
ABSTRACT 

There is currently no explanation of why the corona has the temperature 
and density it has. We present a model which explains how the dynamics of 
magnetic reconnection regulates the conditions in the corona. A bifurcation in 
magnetic reconnection at a critical state enforces an upper bound on the coronal 
temperature for a given density. We present observational evidence from 107 
flares in 37 sun-like stars that stellar coronae are near this critical state. The 
model may be important to self-organized criticality models of the solar corona. 

Subject headings: Sun: flares — Sun: corona — Sun: activity — stars: flare - 
stars: coronae — stars: activity 



INTRODUCTION 



Dynamics in the solar corona takes on a wide range of forms. On one hand, the corona 
is the setting f or the most violent erup tions in the solar system: solar flares and coronal 



mass ejections (lAschwanden et al.ll200ll ). On the other, coronal heating makes the coron a 



almost a thousand t imes hotter than the photosphere, even in the quiet sun (lKlimchukl|2006l ) . 



Parker! (119831 . Il988l ) unified these two phenomena by proposing that micro- and nano-flares, 



studies showing that solar flares exhibit 


power law statistics ( 


Lin et a 


. 1984 




Dennis 


1985; 


Crosbv et al. 1993: Feldman et al. 1997: 


Wheatland 2000: 


Nita et al. 


2002; 


Paczuski et al. 


2005 


) over a wide range of scales for many quantities. [See 


Charbonneau et al 
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2001 


for a 
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review.] In addition, stellar flares ha ve similar light curves to solar flares (Gershbera 12 005) 
and a lso exhibit power law statistics (ICollura et al.lll988l ; IShakhovskayal Il989l ; lAudard et al. 
20001 ). suggesting that the physics of the solar corona is generic to sun-like stars. 



Coronal dynamics rem a ins an active researc h area (jHudsonlll99ll ; iGeorgoulis et al.lll998 



Shibata fc Yokoyamal 120021 ; iHughes et al.l 120031 ). Details of the eruption process including 
how magnetic energy is stored, how eruptions onset, and how the stored energy is converted 
to other forms are still open questions. In addition, while micro- and nano-flares are believed 
to be a major contributor to coronal heating, the authors know of no theory which explains 
why the coronal temperature and density have the values they have, as opposed to larger or 
smaller values. 

In this paper, we pro pose that the condition of the corona is regulated by magnetic re- 
connection (jCassakll2006l ). a dynamical process which converts magnetic energy into kinetic 
energy and heat and energizes particles. Magnetic energy is stored during collisional (slow) 
reconnection, which has been shown to drive the coronal plasma toward lower collisionality 
(ICassak et al.l 120061 ) . If the plasma becomes marginally collision less, a bifurcation in the 
underlying dynamics of reconnection occurs (ICassak et al.ll2007bl ). This bifurcation, which 
occurs when two length scales 5sp and pi (to be defined below) are comparable, catastrophi- 
cally initiates fast (Hall) reconnection, releasing the stored energy in the form of an eruption. 
The condition of marginal collisionality, therefore, sets an upper bound on how hot the coro- 
nal plasma can be for a given density. The continual driving toward lower collisionality of 
the pre-flare corona by slow reconnection enforces the self-organization of the corona to a 
state of marginal collisionality where 5sp ~ Pi- We present finer details of this process below. 
Then, we perform the first observational test of this model using a large sample of data from 
stellar flares on sun-like stars. We find that 8$p and pi are comparable for every event in the 
sample, indicating that stellar coronae do self-organize into a marginally collisional state. 



2. THEORY 



Magnetic reconnection d epends strongly on the col lisionality of the plasma. Collisional 
(Sweet-Parker) reconnection (Sweet 1958 ; Parker 1957h is exceed ingly slow. The thickness 
5sp of the Sweet-Parker diffusion region is given by (jParkerl 119571 ) 



6, 



SP 



7/cr 
4vrc A " 



'SP, 



where ca = B / (Anmin) 1 ! 2 is the Alfven speed, B is the strength of the reconnecting magnetic 
field, rrii is the ion mass, n is the density, rj is the resistivity, and L$p is the length of the 
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Sweet-Parker diffusion region in the outflow direction. The normalized reconnection rate 
Vin/cA ~ 8sp/Lsp is 10~ 7 for coronal parameters, where v in is t he inflow speed. Collisionless 



(Hall ) reconnection has a reconnection rate of the order of 0.1 (jShay et al.lll999t iBirn et al. 



20011 ) . six orders of magnitude faster than Sweet-Parker reconnection for coronal parameters. 



Recent studies (ICassak et al.l 12005k l2007af ) showed that the transition from collisional 
to collisionless reconnection is catastrophic, occurring when 5$p becomes smaller than the 
ion gyroradius pj. At this scale, magnetohydrodynami cs (MHD) breaks down and the Ha ll 



effect (absent in MHD) allows reconnection to be fast (IBirn et a 



2001 



Ro gers et a. 



Cassak et al 



200ll ). 



2007bl ) which 



The transition can be described as a bifurcation (jCassak 20061 ; 
takes stable equilibria out of existence as a control parameter (Ssp/pi ) varies. The releva nt 
gyroradius pi for anti-parallel reconnection is the ion inertial scale di (ICassak et al.l 120051 ) . 



d — — - — 
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where Q C i is the ion cyclotron frequency, u V i is the ion plasma frequency, and e is the ion 
charge. For reconnection with a (guide) magnetic field along the current sheet, th e relevant 
gyroradius becomes p s = c s /Q ci , where c s is the sound speed (ICassak et al.ll2007al ). 



To see how magnetic reconnection self-organizes the corona, consider an active region. 
Before an eruption, the plasma cannot be collisionless: if it were, the stored magnetic energy 
would be rapidly released by Hall reconnection. Therefore, the pre-flare active region must be 
collisional. Since (collisional) Sweet-Parker reconnection is exceedingly slow, magnetic energy 
can be stored. [By "collisional", we mean Ssp > di, which, using eqs. (CQ) and <$fy, is equivalent 
to Vi e > ca/Lsp, where Ui e = i]ne 2 /mi is the ion-electron collision frequency. Therefore, 
reconnection is collisional when the ion transit time al ong the Sweet-Par ker diffusion region 
is longer than the ion-electron collision time. See also lUzdenskyi (j2007bl ).] 



When Sweet-Parker reconnection begins in the corona, e.g., as a result of two coronal 
flux tubes coming together, the reconnecting magnetic field B is initially much weaker than 
the strong asymptotic magnetic field in the core of the flux tube, i.e., the reconnection is 
embedded within a wider current sheet. From equation (JTj), the thickness of the dif fusion 
region will be relatively wide, so that 5sp 3> Pi- It was shown (ICassak et al.l I20061 ) that 
embedded Sweet-Parker reconnection spontaneously self-drives the current sheet to thin- 
ner scales, even without external forcing. This is because the reconnection inflow convects 
stronger magnetic fields into the diffusion region, which causes 5$p to decrease [see eq. ([I])]. 
Thus, the reconnection process itself self-drives the system toward lower collisionality. 



If the asymptotic field is strong enough so that 5$p ~ p«, i.e., the system becomes 
marginally collisionless, then a bifurcation causes Hall reconnection to begin, eruptively 
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releasing the stored energy. [We note in passing that if B in an active region is not strong 
enough to ever satisfy 5sp ~ d% for a given density and temperature then no eruption occurs, 
potentially providing an observational constraint on which active regions erupt and which 
do not.] After the eruption, the corona returns to a collisional state, and the process begins 
again. The continual self-driving of the corona toward lower collisionality keeps coronal 
parameters near the critical condition where the bifurcation occurs (5 s p ~ Pi). 

We propose that this process regulates the temperature of the corona. If the temperature 
T of the corona is larger than the critical value, the Spitzer resistivity 77 is smaller (since 
77 oc T -3 / 2 ). From equation (TjQ), a smaller 77 allows a smaller B to initiate an eruption. 
As such, less magnetic energy is stored and released, and the corona cools. Alternately, 
if the corona is cooler than the critical value, a larger magnetic field is required to set 
off an eruption. More magnetic energy is stored and released, i ncreasing t h e temperatu re. 



Either way, the temperature is driven back to the critical value. lUzdenskyl ( 12006 



2007. 



a.b) 



independently proposed a similar model based on the density, which we discuss in Sec. HJ 



OBSERVATIONAL DATA AND RESULTS 



Observational verification of this model entails confirmi ng; that 5sp and pj are com- 



parab le at fast reconnection onset. Laboratory experiments ( IRen et al.ll2005l ; lEgedal et al. 



20071 ) are consistent with this condition, but direct observations are impossible because the 
length scales are not resolvable. Indirect verification is possible by estimating 5sp and pi 
for coronal parameters. For a solar a ctive region, one finds both length sca l es to be a few 
mete r s, as has been noted prev i ously (IPriest fc Forbes! |2000| ; lUzdenskyl 120031 ; iBhattacharjee 
20041 ; ICassak et al.l 120051 . 120061 ; lUzdenskyl l2007bl ). An important question is whether this 
agreement is indicative of a general mechanism or is just a coincidence for solar parameters. 



We use data from a recent study (IMullan et al.l 120061 ) which analyzed 134 eruptive 
flares from 44 stars of spectral type F, G, K, and M, using the Deep Survey/ Spectrometer 
Instrument on the Extreme Ultraviolet Explorer (EUVE) satellite. The data come from a 
single instrument on a single sat ellite, so ther e are no spurious variations due to different 
instrumental characteristics. See IMullan et al.l (120061 ) for a thorough discussion of the data. 



The e-folding time for the flare signal to decay and the emission measure EM were 
extracted from flare light curves. Because of interruptions in EUVE data during some flares, 
Mullan et al.l (120061 ) presented the Td values for each flare as a unique value, a range of values 
or an upper bound. We retain events having a particular value or a range (using the average) 
and omit events given as bounds, leaving 107 events from 37 stars. 
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Post-flare parameters (the temperature T, density n, minimum magnetic field B m j n , and 



length L and cross sectional area A of coronal l oops) were derived (IMullan et al.ll2006l ) from 
id and EM using an approach due to iHaischl (119831 ) , which assumes that i s comparable 
to the radiative and conductive cooling times. Taking A ~ (L/10) 2 , one finds (jHaischlll983l ) 



0.25^-0.25 



T(K) = a T {EM) 
n(cm" 3 ) = a n {EM)°- 12 % 1A25 
L(cm) = a L (EM)°- 25 T°- 75 



(3) 



where a T = 4 x 10 5 . 



dir. 



10 9 , and ol = 5 x 10 6 are constants (in cgs units). A lower 



bound for the magnetic field B min is estimated by requiring the magnetic pressure B 2 /8ir 
to be at least as large as the gas pressure 2nfc R r, where kp is Bo ltzmann's constant, to 
maintain a coronal loop. As a test of the model, IMullan et al.l (120061 ) surveyed the literature 
for independent measurements of T,n,L, and B min for the stars in their study, finding that 
178 of 212 measurements were consistent with the Haisch model. This justifies treating the 
derived parameters as valid independent of the Haisch model. 

We first verify that the Haisch model gives reasonable result s for solar parameters . 



Eruptive solar flares have ~ 10 4 5 sec and EM ~ 10 49 50 cm 3 (jPriest Forbes! 120001 ) . 



Using these values, the Haisch model predicts post-flare parameters of T ~ (4 — 13) MK, 
n ~ (0.3 - 5.6) x 10 10 cm" 3 , L ~ lO 10 " 11 cm, a nd B min ~ (10- 70) G . Compact solar flares 



have Td ~ 10 3 sec and EM 



10 



47-49 



cm 3 (jPriest fc Forbesl l2000l ). Using these values 



L ~ (0.5 - 1.6) x 
,„,„ are consistent with 



the Haisch model gives T ~ (4 — 13) MK, n ~ (3 — 5) x 10" cm 
10 9 cm, and B m i n ~ (90 — 220) G. These ranges of T 1 n ) L 1 and B 
inde pendent empirical va l ues obtained from images a nd X-ray data for flaring loops in the 
sun (IFeldman et al.l 1 19951 ; IShibata fc Yokoyamal 120021 ) . 



To calculate 5sp and di from the Haisch model, we use B m i n for the upstream magnetic 
field and n for the density. Sweet-Parker current sheets extend to system scales (IBiskamp 
19861 ). so Lsp is on the order of the coronal loop radius A 1 / 2 ~ L /10, consistent with th e 
Haisch model. Finally, we use T to calculate the Spitzer resistivity (jSpitzer fe Harmlll953l ) 



V 



16v^Fe 2 In A 

3m P 



2k B T 



3/2 



(4) 



where m e is the electron mass and InA = ln[(3/2e 3 )(fc 3 3 T 3 /7m) 1 / 2 ] is the Coulomb logarithm. 
Use of this formula is justified because the electron mean free path (A m f Pje ~ v t h,e/ v e i ~ 25 
km for solar conditions, where Vth,e is the electron thermal speed and f e , is the electron-ion 
collision frequency) is small compared to length scales in the outflow direction (Lsp ~ 10 4 
km) and along the current sheet (L ~ 10 5 km). 
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The result of comparing 5,sp to dj using the stellar flare data is plotted in Fig. [TJ 
Representative solar values based on = 10 4,5 sec and EM = 10 49 ' 5 cm~ 3 for eruptive flares 
(Ssp ~ 110 cm and di ~ 200 cm) and r<2 = 10 3 sec and EM = 10 48 cm -3 for compact flares 
{fisp ~ 44 cm and di ~ 35 cm) are plotted as the asterisk and plus, respectively. A dashed 
line with slope of unity is plotted. The agreement is extremely good. A least squares analysis 
gives a best fit slope of 0.98 ± 0.02 with a correlation coefficient of 0.981. 

It is encouraging that the slope of the line in Fig. [T] is consistent with unity. However, 
there are ambiguities in the data analysis. For example, we used di as the critical length 
scale, whereas p s is more applicable to the corona (but more difficult to estimate). These 

1/2 

scales differ by a factor of /? fot , where f3 to t is the ratio of gas pressure to total magnetic 
pressure. If (3 to t ~ 0.1 in the corona, this introduces a factor of a few. The present analysis 
does not intend to distinguish between the two gyroradii; rather, the results demonstrate 
that 5,sp is within a factor of a few of the critical length scale pi in active stellar coronae. 

A caveat of the result in Fig. [1] pertains to how the parameters are derived in the 
Haisch model. Using equations ([TJ, d2J), and fll]), using L S p = L/10, and eliminating B by 
defining the ratio of gas pressure to magnetic pressure in the reconnecting magnetic field 
as I3 rec = 2nk B T/(B 2 /87r), we find (5 S p/di) 2 ~ (e 4 In A/15k 2 B )(27rm e (3 rec /m i ) 1 / 2 (nL/T 2 ). 
Treating (3 rec as a fixed parameter and eliminating T, n, and L using equation (j3J) gives 

where a = (e 4 /15fc|)(27rm e /m i ) 1 / 2 = 1.09 x 1(T 8 cm K is a constant. The slow dependence 
on Td/EM significantly suppresses scatter in the observational data when evaluating Ssp/di. 
However, the magnitude of Ssp/di is unconstrained by the Haisch model, so the slope of the 
line in Fig. [T] being of order unity is significant. Furthermore, since the data obtained using 
the Hai sch model agrees wi th independent determinations of the same quantities from other 



studies ( IMullan et al.ll2006l ). it is reasonable to assert that data obtained independently from 



the Haisch model would fall close to the same line. 

We can avoid suppression of the scatter in the data by solving equation (jSJ) for r d 
and taking a logarithm of both sides. This yields log(rd) = log(EM) + C, where C = 
16\og(8sp/di) — 4 log p rec — 47 using a value of In A ~ 22, which is representative of the 
stellar data in our study. If 5sp ~ di, this predicts a linear relationship between log(rd) and 
log(EM), with C being the y-intercept. The stellar data are plotted in Fig. [21 The gray 



boxe s show the range of values for eruptive and compact flares on the sun (IPriest fc Forbes 



20001 ). Assuming 5sp ~ di and taking (3 rec to be of order unity, the predicted line is plotted. 
While the data do not fall on a line, the line predicted by the hypothesis that Ssp ~ di does 
pass through the data. To see why this is significant, note that if 5$p was, say, 100<ij (at 
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1 — 10m, still a very small length scale compared to coronal loop radii), then C would be —15 
instead of —47 and the line in Fig. [2] would lie 32 units higher, orders of magnitude removed 
from the data. The hypothesis that 5$p ~ d% brings significant ordering to the data. 

We note that the theory predicts 5sp ~ d{ at flare onset, while the Haisch model refers 
to post-flare conditi ons. Following an eru ption on the sun, temperatures typically increase 
by a factor of a few (IFeldman et al.lll995l ). while the density incre ases due to chromos pheric 
evaporation by at least a fact or of 10 [compare pre-flare data (jSchmelz et al.l 119941 ) with 
post-flare data (IDoschekl Il990l )]. From the relation above equation ([5]), Ssp/di oc n l l 2 /T. 
For an increase in n by a factor of 10 and T by a few, 5sp/di does not change appreciably. 
Assuming this is true for other stars, if Ssp ~ di after a flare, it is also true before a flare. 



DISCUSSION 



The data analyzed in this paper pertain to flares in sun-like stars, but the underlying 
dynamics of reconnection is general. Our model applies equally well to micro- and nano- 
flares in the quiet corona. Using values for the quiet sun of T ~ 1 MK, n ~ 10 9 cm -3 , B ~ 5 
G, and L ~ 10 10 cm, we find 5sp ~ 770 cm and d{ ~ 720 cm, in agreement with the model. 

The present result may have important implications for self-organized criticality (SOC) 
models of the solar corona. SOC occurs in driven, dissipative syste ms when the sys tem is 
driven to a critical state where it undergoes a major reconfiguration (IBak et al.lll987l ). SOC 



leads to power law statistics, which encouraged iLu fc Hamilton 



unde r goes SOC. Subsequent stud i es of SOC in the c orona exist (ILu et al.lll993l ; IVlahos et al. 



19911) to propose the corona 



19951 ; lLongcope fc Noonanl I200CH : llsliker et al.ll200ll ). but a firm physical foundation of the 



mechanism for self-driving and the physical condition setting the critical state is often t raded 
for the ease of performing cellular automaton simulations [see lCharbonneau et al.l (120011 ) for a 
review]. The present result provides a physical mechanism for self-driving (embedded Sweet- 
Parker reconnection) and the critical state (marginal collisionality) , which may provide an 
avenue for developing quantitative predictions of SOC to compare with coronal observations. 

An alternate mechanism ( Uzdensky 2006 . 2007al lbl) for heating the solar corona uses a 
change in density to achieve self-regulation. After an eruption, chromospheric evaporation 
increases the coronal density, decreasing the ion gyroradius [eq. (T5])] and making subsequent 
eruptions more difficult. The extent to which Uzdensky 's and our mechanisms regulate 
coronal heating is an open question. 



The present model assumes that Sweet-Parker scaling is appropriate for thin current 
sheets of large extent. Long current sheets are known to fragment due to secondary instabil- 
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ities, but the effect of this on the reconnection rate is unknown. Verification of the present 
model would entail testing whether Sweet-Parker reconnecti o n in ex tended current sheets 
remains much slower than Hall reconnection. [See lUzdenskyl (j2007bl ) for further discussion 
of this point as well as other future research directions.] 

The authors thank J. F. Drake, A. Klimas, E. Ott, S. Owocki, P. So and D. Uzdensky 
for helpful conversations. This work was supported in part by the Delaware Space Grant. 
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Fig. 1. — Sweet-Parker current layer thickness 5$p vs. ion inertial length di = c/u pi for the 
stars in the sample. The dashed line displays their predicted equality. Boxes denote F, G, 
and K stars; diamonds denote M dwarfs. The asterisk at di ~ 200 cm and the plus at di ~ 40 
cm denote v alues based on ayerage and EM values for eruptive and compact solar flares, 
respectively (IPriest &: Forbes! |2000| ). 
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Fig. 2. — Decay time Td vs. emission measure EM for the stars in the sample. The dashed 
line shows the prediction of the theory. Boxes denote F, G, and K stars; diamo nds denote 
M dwarfs. Ranges for eruptive and compact solar flares (IPriest fc Forbes! l2000l ) are shown 
by the gray boxes. 



